Interface dynamics at the depinning transition 
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We study the local scaling properties of driven interfaces in disordered media modeled by the 
Edwards- Wilkinson equation with quenched noise. We find that, due to the super-rough character 
of the interface close to the depinning transition, the local width exhibits an anomalous temporal 
scaling. It does not saturate at a characteristic time ts {I) ~ as expected, but sufi^ers an anomalous 
temporal crossover to a different time regime t^* , where /?« ~ 0.21. This is associated with the 
existence of a local roughness exponent aioc — 1 that differs from the global one a ~ 1.2. The 
relevance of the typical size of pinned sites regions near the critical point is investigated and the 
definition of the critical depinning is dicussed. 
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I. INTRODUCTION 

Rough surfaces and interfaces appear in many situa- 
tions in nature. Deposition, erosion, fluid-fluid displace- 
ment in porous media or flre front motion, are important 
examples [Q,^ in which an interface kinetically becomes 
rough. A rough interface in dimension d + \ (d is the 
substrate dimension), which is described by its height 
h{x,t) at position x and time t, has been treated as 
a self-afhne fractal. In the case of interfaces driven by 
spatio-temporal noise the temporal evolution of /i(x, t) is 
given by the stochastic equations of Edwards- Wilkinson 
(EW) Ig] and Kardar-Parisi-Zhang Q], which have been 
widely studied. The effect of quenched disorder {i.e. 
a noise frozen in time) on the interface dynamics is a 
more difficult problem and up to now there is not a com- 
plete and consistent picture. An important feature of the 
quenched disorder is that the interface can be dominated 
by pinning forces, which can slow down the motion of 
the interface in large regions. Pinning phenomena have a 
considerable importance in problems like the immiscible 
displacement of fluids in porous media or the motion of 
domain walls in magnetic systems. 

During the last years much effort has been paid to this 
subject, but theory, simulations and experiments are not 
yet in good agreement (see |2j for recent reviews). In 
this paper we wish to show that the reason for this dis- 
agreement is certainly twofold. On the one hand, growth 
processes do not always generate truly self-affine inter- 
faces and different scaling behaviours may appear 
at short and large length scales. On the other hand, the 
quenched disorder introduces a new length scale, related 
to the size of pinned regions that becomes relevant and, 
as we will see later, changes the usual scaling laws. 

The simplest stochastic differential equation that de- 
scribes the evolution of the interface height in the pres- 
ence of quenched randomness is the so-called quenched 
Edwards- Wilkinson (QEW) equation: 



which is similar to the EW equation but with a quenched 
disorder term, ?7(x, h). The external driving force F con- 
trols completely the dynamics of the front. If F is larger 
than a critical force F^,, the interface moves with a finite 
velocity. However, the interface remains pinned by the 
disorder for F < Fc. The critical point F = Fc is known 
as depinning transition. Above the depinning transition, 
F > Fc, a. characteristic length, ^, appears. ^ represents 
the typical size of the pinned regions of the interface and 
close to the transition scales a.s ^ {F — Fc)~'^ . As we 
will discuss later, ^ is very important to understand the 
scaling of the interface near the depinning transition. 

The global interface width, which characterizes the 
scaling of the advancing front, is defined by 



a(L,t) = {((Mx,i)-(Mx,i)))2)i/2| 



(2) 



where brackets are averages over the whole system and 
{...} over realizations of the disorder. In the usual case, 
when the only relevant length is the system size, the 
width scales with time for the early times regime and 
saturates at long times: 



aiL,t) 



if t < is 
if i > is 



(3) 



dh 

— = DV^h + F + rjix,h), 



(1) 



where a is the roughness exponent and /3 is the time ex- 
ponent. This scaling picture is known as Family- Vicsek 
scaling ansatz QH]. The saturation time tg depends on 
the system size L because it is the time in which the hor- 
izontal correlation length lc{t) ~ t^^^ reaches the system 
size. The exponent z is called dynamic exponent and the 
scaling relation a = z/3 is fulfilled. 

An alternative way of determining the roughness ex- 
ponent a is frecuently used in simulations because of the 
long wait needed to get saturation in large systems, in 
which simulations are done. The width is calculated over 
a window of length I <C L and averaged over many pieces 
with the same length along the interface. In other words, 
we can study the scaling of the local width 
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the local and global widths. To start we have to discretize 
the QEW equation in the horizontal direction 



The local width scales in a similar manner as the global 
width and has a short time regime a{l, t) ^ . For small 
length scales (or equivalently at long times), I <C lc{t), 
the local width is believed to be independent of time and 
to scale as 



a{l,t) ~ 1° 



(5) 



This method is very useful when one has no possibility of 
producing interfaces in different system sizes as occurs in 
real experiments. Most theoretical growth models gener- 
ate self-affine interfaces and a ~ aioc 

Much work has been carried out to obtain the criti- 
cal exponents at the depinning transition. The rough- 
ness exponent has been the most controversial. Experi- 
ments performed in d = 1 ||^, which are believed to be 
described by Eq.(|l|), gave roughness exponents around 
0.75-0.9 in disagreement with a renormalization group 
prediction a = 1 [^. A scaling theory |^ gave an effec- 
tive roughness exponent 0.8 in good agreement with the 
model reported in Ref . ||l0|] . A number of numerical mod- 
els for interface pinning by quenched disorder have been 
investigated. Roux and Hansen found two different 
exponents: ajoc = 0.86 ± 0.03 when the local width was 
used, and a ~ 1.16 — 1.20 looking at the scaling of the 
global width as a function of the system size. Results 
reported by Jensen |l2j, aioc — 0.9 and a ~ 1.15 are sim- 
ilar. Makse and Amaral jl^] obtained aioc = 0.92 ± 0.04 
or a ~ 1.23 ± 0.04 depending on the method that they 
used. Also other numerical determinations of aioc gave 
values around 0.8 |l|-|l| and a ~ 1.2 

It has been demonstrated that the reason for the exis- 
tence of two different roughness exponents is the super- 
rough character of the interface close to the depinning 
transition. In Ref. |19| was already shown that the local 
width is bounded by (j{l,t) < I and for that, it is just 
technically imposible to measure a roughness exponent 
larger than one using the local width. On the contrary 
the global width does not have this restriction so that it 
gives the true roughness exponent. However, we would 
like to note here that there also exist other effects [|[^] 
that may lead to the same anomalous scaling in surfaces 
even with a < 1. Our aim in this paper is to show that 
the local width suffers a temporal crossover, instead of 
saturating, associated with a different interface scaling 
at short, I ^ L, and long, I ~ L, length scales. The 
scaling properties of the interface cannot be correctly de- 
scribed by the usual dynamic scaling hypothesis. We will 
see that a super-rough dynamic scaling gives the correct 
scaling of the front. 



II. NUMERICAL RESULTS 

We begin by performing a numerical integration of the 
QEW equation in = 1 to study the scaling behaviour of 



h,{t + A) = h,{t) + A{D{h,+i{t) + h, 
+ AiF + r,it,Mt)])), 



,(t)-2h,{t)) (6) 



where the index i = 1, • • • , L, being L the size of the sys- 
tem, [hi (t)] represents the integer part of hi (t) . ho (t) = 
hi{t) and — h^it) are the boundary conditions. 

This is the usual discretization scheme used by many 
authors jl^Jl^,|l^. The random field r]{i,j) is Gaussian 
and its correlation (ri{i, j)r]{i' , j')) = gSi^i'Sj,ji. Hence 
the noise is uncorrelated from a cell to another one and 
stands for the time of a run. The time step A (we used 
A — 0.01 in all our simulations) and the size L have to be 
chosen in such a way that the QEW equation is obtained 
in the continuous limit (A — > 0, L — » oo). So, the follow- 
ing conditions must be verified: AD <C 1, AF <C 1, and 
Ag <C 1. One has to check that changing these parame- 
ters does not affect the final numerical result. Following 
Kessler et. al. we distributed disorder sites randomly 
with a prefixed density p over the network. So, the dis- 
order is nonzero with probability p and on these sites 
disorder is distributed according to a Gaussian distribu- 
tion. The scaling exponents do not depend on the actual 
value of p but the crossover effects will be clearer for small 
values of p. As we were looking for crossover regimes, we 
used p = 0.1 in all our numerical simulations. We worked 
with networks of size L = 1000, and the results were av- 
eraged over 15-50 realizations of disorder. As expected 
we found a critical driving force Fc that separates the 
pinning and depinning phases. The value of Fc depends 
on the actual values of the parameters of simulation (p, 
D, L). For p = 0.1, A = 0.01, D = 5, and L = 1000 we 
estimated Fc ~ 0.13. 



A. Scaling far from the depinning transition 

Our numerical results indicate that the global width 
has two different time regimes. Close to the transition, 
F ~ Fc, the width behaves as a{L,t) ^ t^ , where 
f3 ~ 3/4 is the time exponent near the critical point. 
Far from the pinning transition, F ^ Fc, (j{L, t) ~ t^'^^ 
and the time exponent Pew — 1/4 corresponds to the 
time exponent of the EW equation. As we will see later, 
the values of the remaining exponents a and z indicate 
that the interface dynamics clearly belong to the EW 
universality class in the limit of large driving forces. 

In general, the temporal behaviour of the global width 
is well understood near and far from the depinni ng transi- 
tion and it has been analyzed by many authors ||l^,|l^,|l3| . 
On the contrary, we focus here on the scaling behaviour 
of the local width. Firstly, as far as time evolution is 
concerned, in the strong pushing regime the local width 
scales with time for early times and saturates at time ta{l) 
as corresponds to the Family- Vicsek scaling picture (see 
Fig.l). The value of the width at saturation can be used 
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to determine the roughness exponent. We found that the 
saturated local width, asat{l), scales as asat{l) ~ ^"'^ (see 
Fig.l inset). Thus, for large driving forces the exponents 
(/3 ~ 1/4 and aioc — 1/2) of the EW universality class 
are obtained, as expected ]lO|,p|,p^. 



B. Local scaling at the depinning transition 

At the critical value, F = Fc, the complete scaling de- 
scription of the local width is much less trivial than in the 
strong pushing limit discussed above. In Fig. 2 we display 
the time evolution of the local width a(l,t) for different 
window sizes I. After an early time regime, the local 
width crosses over to t^' , where the asymptotic time ex- 
ponent is /3* ~ 0.25. This scaling region becomes shorter 
as the length scale is larger, I — > L. Moreover, there is 
no saturation of local widths at times ts{l) ^ P, con- 
trary to what is generally thought saturation of the local 
width occurs only when the whole system saturates {i.e. 
ts{l) ^ for any I). As we did in Fig.l for the strong 
pushing limit, we can use the saturated values of the local 
width in Fig. 2 to determine the local roughness exponent 
(see Fig. 2 inset), (Jsat{l) ~ f''^^ . From Fig. 2 it is clear 
that the dynamic behaviour of the interface near the de- 
pinning transition cannot be completely described with 
a Family- Vicsek scaling ansatz. According to this scal- 
ing picture the local width a{l,t) should saturate when 
the horizontal correlation length reaches the size of the 
sample, lc{ts) ^ I, as occurs in the strong pushing limit 
(see Fig.l). However, the local width displayed in Fig. 2 
scales with time for length scales I <C lc{i) and a new 
time exponent (3^ has to be introduced. 

This anomalous time regime of (j{l,t) has been ob- 
served in other growth processes and termed the problem 
of anomalous kinetic roughening in the literature . In 
Ref. 1^ we have shown that for a broad class of growth 
processes the scaling for I <^ lc{t) is actually given by 



(7) 



when a local exponent aioc ^ a exists, i.e. when anoma- 
lous roughening occurs. This gives rise to an anomalous 
temporal regime with [3^ — [3 — aioc/z. Only when the in- 
terface is self-affine one has /?* = and the Family- Vicsek 
scaling for the local width is recovered. As demonstrated 
m Ref. D this may occur for growth with a below or 
above one and thus it is not necessarily associated with 
super-roughening. More precisely, super-roughening cor- 
responds to an anomalous scaling like (|^), but aioc — 1 
and so = a — 1/z in this particular case. At true 
saturation, when lc{t) ^ L, the local width scales as 



a{l,t^ L^) - 



that corresponds to the anomalous dependence of the 
local width on the system size already observed in simu- 
lation Ol. 



Next, we are going to show that for the QEW model 
(|l|) the anomalous scaling behaviour of the local width 
at the depinning transition is also given by Eq.(0) with 
ck/oc = 1 in this case. In Figure 3 we plot a collapse of 
the data of Fig. 2. An inspection of Fig. 3 reveals that 
the dynamic scaling form of the local width is consistent 
with 



a{l,t)^tf'f{l/t^/^). 
The scaling function is given by 



^, ^ I const 
/(«) ~ 1 



if u > 1 
if M < 1, 



(8) 



(9) 



where a/oc — 0.92, (3 ~ 0.81 and z ~ 1.53 are the expo- 
nents giving the best data collapse. This scaling behavior 
leads to a local width a{l,t) - t/3-"ioc/z pioc the in- 
termediate regime I <^ t^/^ ^ L that corresponds to the 
scaling proposed in (|^) with [3^, = (3 — aioc I z — 0.21 (to 
be compared with our previous determination ~ 0.25). 
In agreement with other authors [^3 20 1^ , we thus found 
(3 ~ 0.81, z ~ 1.53 and a global roughness exponent 
a = zf3 ^ 1.2 from the data collapse shown in Fig. 3. 

Similar results were obtained when the structure factor 
was used. The structure factor, which is related to the 
Fourier transform of the height-height correlation func- 
tion can be used to obtain the critical exponents 

S{k, t) = (h(k, t)h{~k, t)) - fc-(2"+i)s(/j2t)^ (10) 

where h{k, t) = L^^l"^ ^YliJ^^^^ ^ ^('•)] cxp(ifca:;) and the 
scaling function is 



s(u) 



const, if u ^ 1 

1*2"+! if u < 1, 



(11) 



whenever the interface satisfies a standard Family- Vicsek 
scaling. Also in the case of anomalous scaling due to 
super-roughening, Eq. (Il^) gives the correct scaling of the 
structure factor (see [5brfor details). 

We have determined the structure factor at the depin- 
ning transition by performing simulations of Eq.(|l|) in a 
system of total size L = 128 and a density of pinning 
sites p — 1. For this system size, the critical force was 
found to be Fc ~ 0.067. Figure 4 shows the data col- 
lapse for the exponents a — 1.25 and z — 1.67 (and then 
(3 = a/z ~ 0.75) in agreement with the exponents that 
we found in Figure 3. 



C. Dynamic exponent at depinning 

There is another point that merits some discussion re- 
lated to the complete scaling of the QEW model. From 
Figs. 2 and 3 we have obtained a value for the dynamic 
exponent zRil.5 — 1.6in agreement with several pre- 
vious simulations |^,^ of the QEW equation (and re- 
lated models) in which an exponent z ~ 1.5 was reported. 
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However, in these numerical works z was always calcu- 
lated from the scaling relation z = a//?, as we have done 
in this paper. In principle, in order to correctly deter- 
mine the dynamic exponent z one should be able to use 
its definition in lc{i) ^ t^l^ ^ where the correlation length 
Icii) can be obtained from an appropriate correlation 
function. Let us consider the height correlation function 
V{l,t) = {{{h{x^l,t) - {h)){h{x,t) - {h)))} , which re- 
lates the heights at two positions x and x + I. This cor- 
relation function becomes zero at distances I larger than 
the correlation length, thus Idt) can be determined from 



lc{t) = 



JlT{l,t)dl 
J T{l,t)dl 



,1/z 



(12) 



Note that Eq.(p^ should be the correct way of obtain- 
ing z while we are interested in checking the validity 
of the dynamic scaling ansatz. Direct determination of 
the dynamic exponent from Eq. (|l^ ) was reported in Ref. 

where z' ~ 2 for a wide range of driving forces 
even at the transition (here the prime denotes a direct 
measure using Eq. (p^) ) . Surprisingly, the value of the 
exponent z' ~ 2 obtained in this way differs from the one 
that we have measured by an indirect method z ~ 1.5 
(and consecuently the corresponding global roughness ex- 
ponents a ~ 1.2 and a' ~ 1.5 are different as well). A 
value of the dynamic exponent z' ~ 2 and the time expo- 
nent P ~ 0.75 leads to /3* — (3—1/z' ~ 0.25, in agreement 
with our previous estimation for the anomalous time ex- 
ponent. 

In order to try to explain this puzzeling situation we 
have considered the effect of the characteristic size of 
the pinned regions on the interface local scaling. Besides 
the super-rough character of interfaces generated by the 
QEW equation and the corresponding anomalous scaling 
of the local width just discussed above, there is also an- 
other important effect that conspires to complicate even 
more the QEW problem. As we have already mentioned, 
the quenched disorder in Eq.(|l|) can slow down the mo- 
tion of the interface in large regions of length ^ in which 
the interface temporarily remains pinned. In an infinite 
system, the size of the pinned regions, ^, diverges as 
^ ^ {F ~ Fc)~'^ when approaching the depinning thresh- 
old Fc from above and the exponent i> ~ 1.35 ||T^ (or 
i/ ~ 4/3 in Ref. jl^) can be measured in simulation from 
a direct observation of these pinned regions in the inter- 
face. For a system of finite size L the depinning transition 
occurs at a value of the driving force F — Fd{L) for which 
{Fd — Fc)~'^ ^ L. As a consequence, in a finite system 
the depinning force Fd{L) depends on the system size, 
and close to the depinning transition both lengths ^ and 
L become relevant in the dynamics. Thus, in any proper 
numerical determination of the global roughness expo- 
nent at the depinning threshold one must have in mind 
two important facts. First, to obtain the global exponent 
a' (where a' = z'/3 « 1.5) one has to perform simulations 
in systems with different sizes L and, after saturation for 
t ^ L^, one can fit the power law (t{L, t L^) ~ L" . 



Secondly, since the critical force Fd{L) is different for a 
different size of the system, in order to keep the system at 
the depinning transition one must vary the driving force 
F in such a way that ^ ~ L. These requirements lead to a 
dependence of the saturated global width on the relevant 
lengths ^ and L at the depinning transition as 



a(L,i»L-)^i"'(|)« 



{Fd{L) ~ F,)--" , (13) 



where 9 is an exponent that characterizes the dependence 
of the global width on the correlation length ^. It is worth 
remarking here that only when ^ ^ L the system is cer- 
tainly being maintained at the depinning transition and 
the exponent a' can be measured. On the contrary, nu- 
merical simulations done in systems with different sizes, 
but in which the driving force F is not adjusted to the 
corresponding Fd{L), always yield a different roughness 
exponent a = a' — 6. The exponent a would appear 
for any F somewhat larger than Fd and gives the global 
roughness in the moving phase. 

In Ref. ||l^ an anomalous dependence of the prefac- 
tor of the width on the driving force as {F — Fc)^'^ was 
already found, with (j> ~ 0.44. We believe that the depen- 
dence of the saturated width on the driving force close 
to the depinning transition indicates that the system is 
not just at the transition. Taking into account the mea- 
sured values ~ 1.35 |l|,|l| and cj) ~ 0.44 |l|l we have 
9 — <j)/v ~ 0.33. So, we obtain a global roughness expo- 
nent a — a' — 9 1.17 in agreement with our numerical 
determination {a ~ 1.2, see Fig. 3 and Fig. 4) as well as 
with previous simulations [ pT|p^ , p^ , p^ , p7t . 



III. DISCUSSION 

In summary, we have found that the complete scal- 
ing behaviour of the QEW equation is not trivial. At 
the depinning transition the interface is super-rough and 
this leads to the violation of the usual scaling ansatz. In 
particular, the local width exhibits a new time regime 
with an exponent /?* ~ 0.21. We have shown that this is 
associated with the existence of two different roughness 
exponents describing the scaling of the local and global 
fiuctuations. A new dynamic scaling, which was already 
used in other growth processes, has been successfully ap- 
plied to understand the scaling behaviour of the local 
width. 

Accordingly to this scaling, Eq. (|^) , the measured value 
of the local roughness, close to unity a/oc — 0.92, suggests 
that the anomalous scaling of the local width is due to 
super-roughening, in such a way that a global roughness 
exponent larger than one must exist. A direct determina- 
tion of the time exponent gives /? ~ 0.75 — 0.81 in agree- 
ment with most of the previous works. On the contrary, 
a direct determination of the global roughness exponent 
from the scaling of the saturated global width with the 
system size can be controversial since variation of the to- 
tal size L leads to a change in the depinning threshold 
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Fc- The effect of the characteristic length of the pinned 
sites regions is relevant close to the depinning transition. 
It seems then that it is possible to obtain two different 
values for the global roughness exponent depending on 
the procedure used. The most usual reported value of 
the global roughness exponent a ~ 1.2 (and z ~ 1.5) is 
obtained by us from collapses of the local width as well 
as from the behaviour of the structure factor. However, 
a direct measure of the correlation l eng ht reflects that 
there is a dynamic exponent z' ~ 2 fl^^ (and conse- 
quently a' ~ 1.5). Moreover, it remains unclear for us 
which of the mentioned methods (and the corresponding 
set of critical exponents) is the proper to describe com- 
pletely the scaling behaviour of the interface at depin- 
ning. We believe that it is a very unsatisfactory situation 
the fact that, on the one hand, the depinning threshold 
changes with the system size and, on the other hand, one 
of the exponents of interest, the global roughness expo- 
nent, has to be measured by changing the system size. 
It seems difficult to define a correct roughness exponent 
at critical depinning as the exponent a in the power law 
ct{L, t — > oo) ~ L" while a variation of L yields a shift in 
the critical point. 

From an experimental point of view the local rough- 
ness exponent is very important since it is the exponent 
accessible in experiments, in which the size of the sample 
remains fixed. We suggest that the anomalous scaling (||) 
might be useful for clarifying whether the QEW equation 
certainly describes the motion of interfaces in disordered 
media, since in many experiments it may be possible to 
study the time evolution of the local width and to deter- 
mine the existence of an anomalous exponent 
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FIG. 1. Time evolution of the local width on diflFerent 
length scales I (window size) for a driving force F = 0.4, 
which is far from the dcpinning transition. Local width sat- 
urates at longer times for larger length scales as corresponds 
to a standard scaling. The values of the width in saturation 
are taken and ploted vs. scale I in a log- log plot (see inset). 
The slope of the line (inset) that fits the numerical data is 0.5 
and gives a local roughness aioc — 0.5 
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FIG. 2. Time evolution of the local width at the depinning 
transition, F = Fc ~ 0.13. The global width, Z = L, is also 
displayed. The continuous straight lines have slopes 3/4 and 
1/4 and are shown to guide the eye. A simple comparison 
with Fig.l reveals the anomalous scaling behaviour for this 
case. The local width does not saturate, but crosses over to 
an anomalous time regime, i^*, where /?» ~ 0.25 and satu- 
ration occurs at the same time for the whole system. Inset 
shows the width in true saturation (flat part in the plots) vs. 
the scale I in a log-log plot. The line (inset) fits the data and 
the slope corresponds to aioc — 0.92. 
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FIG. 3. Data collapse, according to (9), of the results dis- 
played in Fig. 2. The exponents /J = 0.81 and z = 1.53 were 
used. The slope of the straight line gives the local roughness 
exponent aioc = 0.92. 
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FIG. 4. Structure factor data collapse for the QEW equa- 
tion at the dcpinning transition in a system of size L — 128. 
Data shown a good collapse for z = 1.67 and a = 1.25. The 
deviation from the scaling at large abcissa values is due to 
discrete lattice effects. 
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